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We update the theoretical predictions for direct T(nS) hadroproduction in the framework of 
NRQCD. We show that the next-to-leading order corrections in as to the color-singlet transition 
significantly raise the differential cross section at high pr and substantially affect the polarization of 
the T. Motivated by the remaining gap between the NLO yield and the cross section measurements 
at the Tevatron, we evaluate the leading part of the a| contributions, namely those coming from 
T(nS) associated with three light partons. The differential color-singlet cross section at a s is 
in substantial agreement with the data, so that there is no evidence for the need of color-octet 
contributions. Furthermore, we find that the polarization of the T(nS) is longitudinal. We also 
present our predictions for T(nS) production at the LHC. 

PACS numbers: 12.38.Bx,14.40.Gx,13.85.Ni 



Even though the discrepancy between experimental 
measurements and theoretical predictions for T(nS) 
hadroproduction is less dramatic than for J/tp and ip' , 
it is still unclear which are the dominant mechanisms 
at work in those processes pQ. The leading order (LO) 
ctg result in the v expansion of non-relativistic QCD 
(NRQCD) j3j neither predicts the yield seen in the 
data [H [5], undershooting it by a factor of ten at pt 
around 20 GeV, nor shows the correct pt behavior. In 
this case it has been shown that the inclusion of higher- v 
corrections, in particular those associated with the color- 
octet mechanism, allows the data to be described. The 
yield and the px shape can be reproduced, at the cost 
of introducing two unknown long-distance-matrix ele- 
ments of NRQCD that can be determined from fits to 
the data [5J. However, the polarization predicted fol- 
lowing this approach disagrees with both the published 
CDF [5] and new preliminary D0 measurements [7] . 

In view of such a puzzling situation, it is worth reex- 
amining in detail the available theoretical predictions. In 
Refs. [HE], the complete set of a A s corrections to color- 
singlet production was calculated for the first time. The 
results confirmed the importance of the new channels 
that open up at this order, which is due to their different 
Pt scaling (see Fig.[T](c,d)). In this Letter we apply those 
results to a phenomenological study of T production at 
the Tevatron, comparing the differential cross section at 
next-to-leading order (NLO) with the data from the CDF 
collaboration. We also argue that some a| contributions 
coming from three-jet configurations, i.e. T+jjj, such as 
those arising from gluon-fragmentation (Fig. [T] (e)) and 
"high-energy enhanced" (Fig. [I](f)) channels are impor- 
tant. We show that by adding such cv| contributions 
to the NLO yield, the experimental measurements arc 



reproduced. We also predict the polarization of the T 
from the angular distribution of the leptons produced in 
the decay. 
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FIG. 1: Representive diagrams contributing to T hadropro- 
duction at orders a% (a), a% (b,c,d), a| (e,f). See discussions 
in the text. 

We first present the NLO predictions for the pr spec- 
trum of the T at leading order in v, following the same 
procedure as in Ref. [S] and adding information on the 
polarization of the quarkonium state. In our numeri- 
cal evaluation of the differential cross section and po- 
larization parameters, we have used: |i?x(is)(0)| 2 = 
6.48 GeV 3 and |i? T ( 3 s)(0)| 2 = 2.47 GeV 3 pQ; /*„ = 
y/(2m b ) 2 +p T ; Br(T(lS) -► = 0.0248 and 

Br(T(35) -> = 0.0218; m b = 4.75 GeV; PDF 

set: CTEQ6M [10]. For comparison, we also plot the LO 
yield, for which we used the PDF set CTEQ6L1. 
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The experimental results for the cross section differ- 
ential in pt for prompt T(15) and T(25), as well as 
the cross section for direct T(3S) and the polarization 
for prompt T(15) were obtained in Run I at y/s =1.8 
TeV 0. For Run II at =1.96 TeV, so far only D0 
have analyzed data for T(15 l ) [IT]. The comparison per- 
formed in |nj shows full agreement with [5], once the 
difference in y/s is taken into account. 

In order to compare our calculation for direct T(IS') 
production with the CDF data, we have multiplied the 
most recent prompt-T(lS') cross-section measurement 



of by the direct fraction F^ t > = 0.5 ± 0.12. This 
fraction was obtained from an older sample 4 , where 
the selection cuts \pr\ > 8 GeV and \r)\ < 0.7 were ap- 
plied. The similar px dependence observed for prompt- 
T(15) and direct-T(35) production, as shown in Fig. 2 
of [5], confirms that this fraction should not depend very 
strongly on px- This justifies using the same fraction for 
the whole px range. The errors from the prompt cross 
section and the direct fraction have been combined in 
quadrature. Similar information on the direct yield of 
T(25) is not yet available. 

Our results for direct production of T(15) and T(35) 
at y/s — 1.8 TeV are shown in Fig. U The curves 
labeled as NNLO* will be discussed later on. In or- 
der to estimate the theoretical uncertainties, we have 
varied the renormalization and factorization scales be- 
tween 0.5 no and 2/j.q, keeping them equal. We have 
combined the resulting uncertainties in quadrature with 
those resulting from variation of the bottom quark mass, 
nif, = 4.75 ± 0.25 GeV. The uncertainty associated with 
the non-perturbative parameter |i? T ( n s)(0)| 2 [BJ[T5], cor- 
responding to an overall normalization, is not included. 
The curves for direct T(25) production can be easily ob- 
tained - at this order in v - by adopting the correspond- 
ing values for the leptonic branching ratio and for the 
non-perturbative matrix clement, |-Rx(2S)(0)| 2 - 

As a first comment, we note that the LO and NLO re- 
sults are very similar at low px, i.e. where the bulk of the 
cross section is found. This implies that the a| correc- 
tions are small in this region and that the perturbative 
expansion of the total cross section is under control. On 
the other hand, the relative importance of the NLO cor- 
rections dramatically increases at larger px, where new 
p^ 6 channels (Fig. ll(c)) start to dominate. This behav- 
ior significantly reduces the discrepancy with the experi- 
mental data that is found at LO. However, the NLO pre- 
diction still drops too fast in comparison with the data. 
Note that the NLO calculation also takes into account 
associated production T + bb (Fig. [l] (d)). This contribu- 
tion, in spite of reproducing the shape of the data very 
well at both low and high pt (asymptotically it behaves 
as p T 4 ), remains negligible up to at least px — 30 GeV. 

Throughout our calculation, the polarization informa- 
tion of the T can be traced through the angular distri- 
butions of the two decay leptons. Defining 9 as the angle 
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FIG. 2: Results for (a) T(1S)+X and (b) T(3S) +X at LO 
(a|), for associated production (as), for full NLO (a|+ Q s) 
and for NNLO* (up to a.%) compared with the direct yield at 
= 1.8 TeV measured by the CDF collaboration 0. The 
theoretical-error bands for LO, NLO and associated produc- 
tion come from combining the uncertainties resulting from the 
choice of /j,/, fi r , mi- The uncertainty of NNLO* includes the 



variation of the cutoff, 
the scales fi r and fif. 
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FIG. 3: Polarization for direct T(nS) (n = 1, 2, 3) production 
at ^3= 1.8 TeV up to the order a% (NNLO*). Most of the 
uncertainties on a for LO, associated production and NLO 
cancel. The uncertainty band of the (NNLO*) result comes 
from the variation of the cutoff s"" n . 
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between the £ + direction in the quarkonium rest frame 
and the quarkonium direction in the laboratory frame, 
the normalized angular distribution /(cos 9) reads 

J(cos0) = - ' (1 + a cos 2 9), (1) 
2(a + 3J 

from which we have extracted the polarization parameter 
a bin by bin in pr- Our results are shown in Fig. [3] along 
with the curves for the LO yield and for the associated 
production alone. Even though the inclusion of NLO 
corrections brings the predictions closer to the prompt 
data, the discrepancy in the normalization makes such a 
comparison not very compelling. 

As already noted, the discrepancy between the data 
and the full NLO (up to a|) result in Fig. [2] increases 
with pt- This might be taken as an indication that one 
or more processes with a more gentle px scaling than the 
otg contributions, which scale at most as p^ e (disregard- 
ing the sub-dominant T + bb channel), arc still missing. 
In this section, we argue that the a| terms provide us 
with the missing contributions. Indeed, the gluon frag- 
mentation channel (Fig. [T] (e)) occurs for the first time 
at order a| and behaves like p^ at large transverse 
momentum. At next-to-next-to-leading order (NNLO) 
"high-energy enhanced" channels (Fig.[l](f)) also appear 
for the first time, whose contributions have previously 
been approximated in the factorization approach |13j . 
Both these contributions, which exhaust the possibilities 
of new kinematical enhancements at higher order, are at 
the Born level and therefore finite at a|. They are a 
subset of pp — ► T + jjj (j stands for u,d, s, c, g) and pro- 
vide us with new mechanisms to produce a high-p^ T 
with a lower kinematic suppression. They can therefore 
be expected to dominate the differential cross section at 
NNLO accuracy in the region of large transverse momen- 
tum. 

We present here an estimate of these dominant a| 
contributions at large px, bypassing the presently out- 
of-reach inclusive calculation of pp — > T + X at NNLO 
accuracy. We consider the whole set of a| processes 
contributing to the color-singlet production of a T as- 
sociated with three light partons. The computation of 
the associated matrix elements is technically possible via 
the approach presented in Ref. [14] . In order to protect 
the integrated cross section from soft and collinear di- 
vergences, we impose a "democratic" invariant-mass cut 
s"j ,n of the order of ml to any pair of light partons so 
that the phase-space integration remains finite. For new 
channels opening up at a|, the dependence on the value 
of s™ ln is expected to be small in the region pr ^> nib, 
since no collinear or soft divergences can appear there. 
For channels whose Born-level contribution occurs at a| 
or otg, the invariant-mass cut results in potentially large 
logarithms of Sjj/s™ m in the differential cross sections. 
However, these logs factorize over lower order amplitudes 
that are suppressed by powers of px compared to the 



dominant a| contributions. The result is that the sensi- 
tivity to the value chosen for the minimal invariant mass 
cut s™ ln is expected to die away as px increases. 

In fact, this argument can be tested by using it to 
estimate the dominant contribution to the production of 
T at NLO accuracy. The differential cross section for 
the real corrections, T + jj production, is displayed 
in Fig. [4] The grey band is obtained by varying the 
invariant-mass cut s™ m between any pairs of light partons 
from 0.5ml to 1m\. The yield becomes insensitive to 
the value of s™ m as pt increases, and it reproduces very 
accurately the differential cross section at NLO accuracy. 
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FIG. 4: Full computation at NLO for T + X (dashed line) vs. 
T(1S) + 2 light partons with a cut on s™ m (grey band). The 
red band shows the cross section for T(15") + 3 light partons 
with the same cut. (see text for details). 

The computation of pp — > T + jjj at tree level is in 
principle straightforward, but technically quite challeng- 
ing: several parton-level subprocesses contribute, each 
involving up to hundreds of Feynman diagrams. We fol- 
low the approach described in Ref. [14] , which allows the 
automatic generation of both the subprocesses and the 
corresponding scattering amplitudes. We find that the 
subprocess gg — ► T + ggg dominates, providing 50% of 
the whole yield. The uncertainty associated with the 
choice of the cut s™ ln is substantial although, as can be 
seen in Fig. |4j it decreases at large pt where these con- 
tributions are dominant. 

The differential cross-sections for T(15) and T(3S) are 
shown in Fig. [2] The red band (referred to as NNLO*) 
corresponds to the sum of the NLO yield and the T + jjj 
contributions. The contribution of T production with 
three light partons fills the gap between the data and the 
NLO calculation. The a| contribution is very sensitive 
to the choice of the renormalization scale, u, r . This is 
expected: for moderate values of the pt, the missing vir- 
tual part might be important, whereas at large pt, the 
yield is dominated by Born-level a|-channels from which 
we expect a large dependence on /i r . Even though the 
uncertainty on the normalization is rather large, the pre- 
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diction of the pt shape is quite robust and agrees well 
with the behavior found in the data. At leading order in 
v, results for T(2S) can be obtained simply by changing 
\Rr(nS) (0) | 2 an d the branching ratio. The predictions for 
the polarization parameter a are only slightly affected by 
the NNLO* contributions and remain negative at high 
transverse momentum, Fig. [3] Finally, we show our pre- 
dictions for the pt distributions at the LHC (y/s = 14 
TeV) in Fig. [SJ 

Obviously, our analysis cannot be extended to rather 
low pt, where the approximations on which it is based no 
longer hold. In this case, one could improve the predic- 
tions by merging the matrix elements with parton show- 
ers using one of the approaches available in the litera- 
ture |15j , or by performing an analytic resummation |16j . 



NLO (a s ) and the data. 

As a further result, we have predicted the polarization 
parameter a to be negative and only slightly decrease 
with pt, indicating a longitudinally polarized yield. More 
accurate data, in particular for T(3S') and possibly at 
higher pt, are needed to check these predictions. 

In summary, the comparison with the available data 
suggests that there is no need to include contributions 
from color-octet transitions for T(nS) production. For 
the ip case, the situation is less clear since the addition 
of higher-order QCD corrections seems not to suffice to 
describe the data [T71 [TS]. In this case, other mecha- 
nisms may be at work such as those recently discussed in 
Refs. [191 EQ]. 
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FIG. 5: T (IS) + X at yfs = 14 TeV. Same parameters as for 
Fig 2. Results for the T(3S) +X can be obtained by rescaling 
the curves by the ratio of the corresponding wave functions 
values at the origin and branching ratios. 

In conclusion, we have compared the full NLO (a s ) 
cross section for polarized direct hadroproduction of 
T(nS) to the data from the Tevatron. At this order, 
the polarization parameter a is negative and decreases 
with pt- However, the yield is still too low compared to 
the data especially at high pt- 

This led us to consider a subset of the NNLO (a s ) cor- 
rections, namely those corresponding to the production 
of a T associated with three light partons and containing 
p^ 4 contributions. Integrating the amplitudes for such 
processes at a| leads to divergences in some phase-space 
regions that we avoided by imposing a minimal invariant 
mass between the light partons. This was motivated by 
the fact that the neglected contributions are proportional 
to ag amplitudes, which scale at most as p^ 6 , so that we 
expected this approach to give a reliable estimate of the 
full a| result at large pt- We validated our expectations 
by comparing such an approach with the computation 
at oig and by explicitly checking the sensitivity of our 
results to the cutoff choice at high pt- The estimated 
ag-contributions appear to fill the gap between the full 
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